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ABSTRACT 


Low frequency active sonar performance in shallow water is often limited by 
reverberation. Reverberation modeling in shallow water has been difficult due to the 
complexity of the multipath acoustic propagation problem inherent in shallow 
environments. In August 1992, a shallow water, low-frequency reverberation 
measurement was made in the Barents Sea utilizing explosive "signal, underwater sound” 
(SUS) charges as sound sources and a 16-element vertical hydrophone array as the 
receiver. The objectives of this thesis were to analyze the reverberation data from this 
experiment, compare several theones which have been proposed to model reverberation, 
and determine the reverberant characteristics of the region. The three-dimensional 
Hamiltonian Acoustic Ray-tracing Program for the Ocean (HARPO) was used as the 
primary propagation modeling tool. The temporal signal processing consisted of a short- 
time Founer transform spectral estimation method applied to data from a single 
hydrophone. Chapman’s source spectrum model was used. Reverberation models based 
on Lambert’s law and omnidirectional backscattering theory were compared. Lambert’s 
law was found to be more applicable in the Barents Sea. A statistical analysis was 
performed on broadband and narrowband hydrophone data showing that reverberation in 


the Barents Sea possesses Gaussian properties. 
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I. INTRODUCTION 


A. MOTIVATION 

As the world political climate continues to shift, the new naval battlefield has been 
extended into littoral regions (O’Keefe, 1992). Prosecution of a new breed of coastal 
diesel submarines has become a top priority for ensuring U.S. forces maintain battlespace 
dominance in these regions. The development of low frequency active acoustic (LFAA) 
systems by several nations worldwide, including the United States, may provide new, 
more effective methods by which such forces can be combatted. Such systems represent 
a double-edged sword, however, because their potential effectiveness in the detection and 
prosecution of diesel submarines operating in littoral regions implies a threat to the 
tactical secunty of U.S. submarines. Hence, performance predictions for LFAA systems 
have become a priority for the Navy. 

The performance of LFAA systems in shallow water is often limited by 
reverberation from the bottom. Reverberation modeling is a complex task, for it involves 
nearly every aspect of underwater acoustics research. It incorporates environmental 
monitoring, propagation modeling, underwater explosion theory, array analysis, and 
scattering theory into a single field of research. Because the problem is so complex, 


reverberation modeling has been historically difficult due to the lack of adequate 


processing power. Hence, even though the need for accurate predictions has existed for 
many years, only recently has it been possible to adequately model the problem. 

In the past, much research effort has been expended in the study of reverberation 
in deep water. Most of the experiments have involved relatively simple propagation and 
scattering models (McCammon, 1993). They have mainly used point sources and single 
hydrophone receivers to analyze energy which traveled on a single path to a scatterer and 
then returned on that same path. This lessened the computational load significantly, but 
limited the range for which the analysis was accurate. Nonetheless, in deep water, these 
relatively simple models were adequate to produce relatively accurate performance 
predictions of sonar systems. 

In shallow water, due to the multipath acoustic propagation taking place, such 
models were not adequate. Hence, a more accurate, and necessarily more complicated 
modeling technique was needed. This need motivated the inclusion of a low frequency 
reverberation study as part of the Barents Sea Polar Front Experiment (BSPFEX). In this 
experiment, a vertical array and explosive "signal, underwater sound" (SUS) charges were 
utilized to conduct low frequency reverberation measurements. These measurements then 
provided baseline data which have been used in this thesis to test the validity of 
reverberation prediction models. 

The experiment also provided a valuable opportunity to investigate the statistical 
properties of reverberation in the Barents Sea. In a shallow water environment with a 
rigid bottom such as the BSPFEX experiment area, the highly multipath nature of acoustic 


propagation creates a reverberant signal which is a combination of many random 


processes, each of which may or may not be Gaussian. By the central limit theorem, 
however, the sum of these signals arriving at a single hydrophone should possess a 
normal Gaussian distribution. This hypothesis was examined as part of this BSPFEX data 


analysis for both narrowband and broadband signals. 


B. OBJECTIVES AND APPROACH 

The main thrusts of this thesis can be summed up in three objectives. The first was 
to examine two commonly used laws for reverberation prediction as applied to the 
Barents Sea, namely Lambert's law and omnidirectional scattering. The second objective 
was to estimate the environmentally-dependent reverberation model parameter p for the 
Barents Sea and to determine its frequency dependence. The third objective was to 
perform a statistical analysis of the BSPFEX data to determine if the characteristics of 
reverberation in the region did in fact possess Gaussian properties. 

Considering the oceanographic data available from unclassified sources and the 
computing resources on hand, the following methods/models were used for theoretical 
predictions and data analysis in this thesis. Narrowband data analysis was performed 
using short-time Fourier transforms for all power spectral density estimates. Statistical 
analysis was performed on broadband and narrowband data using normalized data 
segments and histograms. Propagation modeling was performed utilizing ray theory as 
implemented in the Hamiltonian Acoustic Ray-tracing Program for the Ocean (HARPO), 
with post-processing performed on the HARPO output to determine the effects of 


spreading and boundary interactions. Geometric effects were accounted for using a ray- 


tube spreading algorithm. Boundary interactions were accounted for using the Rayleigh 
reflection coefficient and Eckart rough surface scattering function. The propagation due 
tO a point source was approximated using four ray launch azimuths ("4 x 2D" estimate). 
Backscattering functions developed by Lambert, Tolstoy and Clay, and McCammon were 
investigated. A SUS source model developed by Chapman (1988) was used. Finally, the 
propagation model, backscattering function, and source model were combined to produce 
a narrowband reverberation envelope prediction for the Barents Sea. By comparing 
theoretical predictions for the Barents Sea with the experimental measurements, the 
validity of each reverberation theory was tested. 

The most significant choice of the the models/methods listed above was that of the 
propagation model, for it had the greatest effect on the level of the predicted signal and 
it determined which models or methods were appropriate for boundary interactions and 
scattering. All routines developed for this thesis were programmed in the MATLAB? 
environment and implemented on an HP 9000 series 735 computer (TAC 3), a general 


purpose UNIX based machine. 


C. OVERVIEW OF THESIS 

The approach described above is summarized in Figure 1.1. Environmental data 
were gathered from unclassified sources, including bathymetric and geologic data from 
recent publications and sound speed field data determined during the BSPFEX. HARPO 
was used to trace rays on several azimuths from the array location in order to estimate 


the three-dimensional propagation effects which determine the transmission paths between 


Our point source (SUS), bottom locations, and a hydrophone. A SUS source spectrum 
model from recent research was used as a second input to the reverberanon prediction. 
The last factor needed was the bottom scattering strength function S,(f,0,0). Scattering 
functions suggested by Lambert, Tolstoy and Clay, and McCammon were considered. 


With these inputs, the reverberation level present at the array was predicted. 
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Figure 1.1: Reverberation level measurement and prediction flowchart. 
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The measurement portion of the problem was addressed in a somewhat more 
straightforward manner. A spectral analysis subroutine was applied to the hydrophone 
data to produce a time series of reverberation level at specific frequencies. Frequencies 
of 50, 100, 200, and 400 Hz were considered in this thesis. These frequencies have no 
special significance (i.e. to an existing LFAA system), other than the fact that they cover 


the array’s 50 to 500 Hz bandwidth in octave steps. The measured and predicted levels 


were compared to determine the validity of the scattering functions used. The statistical 
analysis of the BSPFEX data was performed as previously described and compared with 
a normal Gaussian distribution to verify its statistical properties agreed. 

In Chapter II, details of the reverberation experiment in the BSPFEX are presented. 
The discussion includes a brief discussion of oceanographic and geologic features of the 
Barents Sea, a description of the SUS charges used as sound sources, and the construction 
and operation of the vertical array. This is followed in Chapter III with a discussion of 
reverberation theory, both in general terms and specific theones which have been 
developed. Theories to be presented include Lambert’s Law, Tolstoy and Clay’s low 
frequency scattering function, and McCammon's low frequency, low grazing angle 
backscattering strength model. Chapter IV presents a method for predicting bottom 
reverberation using HARPO, a source model, and a backscattering strength model. It also 
presents the results of the BSPFEX and compares them to the predicted levels. It ends 
with an analysis of the statistical properties of the reverberation measured during the 
experiment. Finally, Chapter V summarizes the work completed and discusses lessons 


learned. Advice regarding future work to refine the present study will be provided. 


IL BARENTS SEA EXPERIMENT 


A. EXPERIMENT OBJECTIVES 

The Barents Sea Polar Front Experiment (BSPFEX) was conducted in August of 
1992. The experiment was a joint effort between the Naval Postgraduate School (NPS), 
Woods Hole Oceanographic Institute (WHO)), and the Science Applications International 
Corporation (SAIC). The principal investigators for the experiment were Professors 
Ching-Sang Chiu, James H. Miller, and Robert Bourke from NPS, Dr. James F. Lynch 
from WHOI, and Dr. Robin Muench from SAIC. The principal engineers for the 
development and deployment of the vertical hydrophone array system were Mr. Keith von 
der Heydt and Mr. John Kemp from WHOI. The objectives of this experiment as 
outlined by the Barents Sea Polar Front Group (1992) were: 

1. Provide a detailed physical description of the polar front. 


2. Enhance the understanding of dynamics of the front, including frontogenesis 
and its influence on regional oceanographic processes. 


3. Assess the ability of acoustic tomography to define frontal and associated 
mesoscale features. 


4. Provide improved documentation of shallow water acoustic propagation in 
this region and the effect of the environment on acoustic Anti-Submarine 
Warfare (ASW) operations. 





LONGITUDE 


Figure 2.1: The Barents Sea. 


The reverberation experiment was designed to meet part of the fourth goal by studying 


low frequency sound propagation and scattering in shallow water. (McLaughlin, 1993) 


B. EXPERIMENT LOCATION 
The experiment was performed in the Barents Sea 100 km east of Bear Island, as 
shown in Figure 2.1. The Barents Sea forms an epicontinental sea (1.3 million km?), 


bounded by the Arctic Ocean to the north, the Svalbard archipelago and the Norwegian- 
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Figure 2.2: Bathymetry in the BSPFEX area. 
Greenland sea to the west, the Fennoscandian shield in the south and Franz Josef Land 
and Novaya Zemlya to the east. Itis characterized by northeast-southwest trending basins 
300 to 500 meters in depth, with an average depth of 230 meters over the entire region. 
It is thus significantly deeper than most of the present day high Arctic shallow shelves 
outside North America, Northern Europe and Northern Asia (which are generally 10 to 
60 meters deep). The greater depth is most likely a response to repeated glaciations in 
the Late Cenozoic, leaving only a thin sediment cover above the Mesozoic and Paleozoic 


bedrock (typically 5 to 15 meters thick) (Norsk Polarinsttutt, 1987). 


The majority of the Barents Sea is covered with a thin sediment layer, though large 
sediment distributions are present in water depths exceeding 300 meters in the western 
part of the major troughs Bjdrndyrenna and Storfjordrenna and exceed 500 meters in 
thickness near the shelf edge. In the experiment area, the sediment thickness varies from 
nearly 0 msec (exposed bedrock) in the northern portion to approximately 50 msec near 
the southern extent of the area. The majority of the sediment is composed of stiff pebbly 
mud, till and/or glaciomarine deposits overrun and reworked by glacial activity (Norsk 


Polarinstitutt, 1987). 
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Figure 2.3: Sound speed profile near the vertical line array. 
Near the experiment area, the oceanography is dominated by the Barents Sea Polar 


Front. It is a stable feature which is caused by circulation in the region (Emblidge, 1991) 
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and which passes directly through the experiment area approximately 3 km north of the 
receiving array. The physical oceanography of the region plays a large role in controlling 
sound propagation. The bottom, as previously described, is rough in some areas and 
consists of a thin sediment overlying hard rock. The bathymety (Figure 2.2) and sound 
speed are highly range dependent. The shallow water and generally downward refracting 
sound velocity profile cause many surface and bottom interactions in the acoustic 
propagation. 

The sound speed field was observed to vary significantly with range. The sound 
speed profile measured at the vertical line array is shown in Figure 2.3 and is a typical 
Barents Sea profile south of the front for the summer. The sound speed gradient is 
basically negative from surface to bottom with a shallow mixed layer region. The 
negative gradient tends to bend rays toward the bottom, enhancing the region’s already 


strong bottom reverberation characteristics and minimizing surface reverberant effects. 


C. EQUIPMENT DESCRIPTION 
The locations for the vertical line array and the SUS charges which were processed 


are shown in Figure 2.4 and are listed in Table 2.1. A description of each follows. 


1. SUS Charges 
A total of 24 SUS charges were expended during the experiment. Eighteen 
of the charges were delivered by a P-3 Orion aircraft which overflew the data recording 
ship on two courses, 242°T and 351°T, dropping 9 SUS charges on each run. The SUS 


charges dropped by the P-3 were U.S. Navy Mark 64-0 exercise communication SUS, 


1] 
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Figure 2.4: Device locations for reverberation experiment. 

which contain a 1.1 lb (0.50 kg) TNT main charge and were set to explode at 60 ft (18 
meters). The charges have a nominal broadband source level of 263 dB re 1pPa (O 1 
meter (U.S. Navy, 1989). All of the air-dropped charges exploded and were recorded, but 
the 60 foot detonation depth placed the charges well within the mixed layer. Although 
they were easily heard by the array, the reverberation created by the charges was minimal 
and the high propagation loss between the detonation points and the array severely 
attenuated the signal. Hence, the recordings of the air-dropped SUS are of little use for 
this particular research. 

The remainigg six charges were thrown by hand off the back of the recording 


ship. One of the six was a dud, and one was not recorded. The ship-dropped SUS 
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TABLE 2.1: DEVICE LOCATIONS 


Device Latitude Longitude Distance to 
Array 


e Line 749 19.1994'N 23° 32.2957'E | 


m WIE Da daas 


SUS 4 74? 18.54'N 23" 3308 E - 1405 meters 


SUS 5 749 18.97'N 239 34 72"E 1287 meters 
SUS 6 749 18.59'N 23? 3420'E 









nb m o o 













charges used were U.S. Navy Mark 82 special purpose SUS, which contain a 1.8 Ib (0.82 
kg) TNT main charge and are set to explode at a depth of 300 ft (91.4 meters). The 
charges have a nominal broadband source level of 274 dB re lpgPa (9 1 meter (U.S. Navy, 
1989). They were dropped at an average range of 1350 meters from the array and 
produced a reverberation signal which is distinguishable above the background noise level 
for approximately 50 seconds. 

Qualitatively speaking, Mk 82 SUS provide an impulsive signal which. 
although non-linear in the immediate vicinity of the source, serves as a very close 
estimate of a point source with an impulse-like pressure signature. A time series of one 
such explosion is shown in Figure 2.5. As noted in the following paragraphs, the signal 
would be even more impulsive in appearance if it had not been for non-linearities in the 
receiving array. A more complete discussion of underwater explosion theory is discussed 


in Chapter III, Reverberation Theory. 
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Mk 82 SUS Level at 1.2 km., Received at Phone 0 of Array 
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Figure 2.5: Received direct blast and reverberation for SUS 3. 


2. Vertical Line Array 
The vertical line array is shown in Figure 2.6. The array consisted of 16 
hydrophones spaced 10 meters apart. The hydrophones had a nominal sensitivity of -160 
dB re V/pPa. The mooring was designed to acoustically decouple the surface buoy from 
the array and hence minimize the effect of surface waves on mooring motion. The 
surface buoy contained an RF Ethernet system which transmitted the digitized data to the 
ship (Von der Heydt, et al., 1992). On the ship the data was stored on tape for later 


analysis. 
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Figure 2.6: Vertical line array and RF Ethernet buoy. (Von der Heydt, et al., 
1992) 


In the tomography experiment, the sign (+/-) of the signal made little 
difference, since the only significant data was the travel time of the signals. In the 
reverberation experiment, it was found that the received signals were consistently 180° 
out of phase with theory due to the details of the data acquisition system. This shows up 
most vividly in the initial blast waveform, which should be initially an overpressure in 
the ground wave, then the waterborne shock wave, again positive-going, followed by 
several positive-going bubble pulses. After the initial arrivals, arrivals from other 
multipaths, including surface reflections (which do undergo a 180° phase shift and hence 
would be initially negative-going) come in but are complicated by the bubble pulses of 


the previous arrivals. The corrected signal is shown in Figure 2.7. In all data analysis, 
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Figure 2.7: Initial received waveform from Mk 82 SUS. 
the polarity of the received signal will be reversed in order to account for this feature. 
An additional undocumented feature of the vertical line array was discovered 
during analysis of the SUS explosion time series. Either the hydrophones on the array 
or the analog-to-digital (A/D) converter in the buoy showed a non-linear response at 
levels exceeding ~160 dB (1 volt), with the signal entirely clipped above 165 dB (1.8 
volts). Figure 2.8 illustrates this clearly during the first half second of the received blast. 
The signal should have been a sharp increase to a peak (at -195 dB), followed by a near 
linear decrease in amplitude (in decibels). This would correspond to an exponential 
decrease in pressure with time. Instead, the signal remains at the maximum phone output 


level until the initial shock wave and first bubble pulses pass. The A/D converter had 
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Mk 82 SUS Level at 1.2 km., Received at Phone 0 of Array 
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Figure 2.8: Reverberation level plot showing clipping of hydrophone output. 


been specifically designed to show a linear response over -5 volts to 5 volts hydrophone 
output. Therefore, the phones themselves were likely the source of the non-linear 
behavior. To be conservative, it is assumed that the phones display a non-linear 
amplitude response over 150 dB re 1 uPa, as is illustrated graphically in Figure 2.9. 
The implication of this on the reverberation experiment is that approximately 
the first second of the SUS explosion data is unusable. This is unfortunate, since any 


time-domain methods of reverberation require a time series as the input signal to the 
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reverberation equation. Because of this, spectrum levels of SUS charges from recent 


research have been used in all data analysis and prediction routines. 
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Figure 2.9: Apparent hydrophone response curve. 
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III. REVERBERATION THEORY 


A. GENERAL REVERBERATION THEORY 

Any discontinuity in the physical properties of a medium tend to intercept and 
reradiate a portion of the acoustic energy incident upon them (Urick, 1983). This 
reradiation of energy is called scattering, and the sum total of scattering from all 
scatterers in the medium to a receiver is called reverberation. These discontinuities can 
be due to any type of inhomogeneity present in the sea. Particles of inhomogeneous 
material (dust, fish, bubbles, etc...), the sea surface, the bottom, and changes in sound 
velocity in the medium away from a perfect isovelocity structure (vertical or horizontal 
effects) all contribute as scatterers. 

Reverberation can be generally categorized into two broad areas: Volume and 
surface. Volume reverberation is due to energy which has undergone scattering with 
Some inhomgeneity which 1s suspended in the water volume and then has propagated back 
to the desired receiver location. Examples of scatterers which contribute to volume 
reverberation include marine life, suspended sediment particles, and bubbles. Surface 
reverberation is due to energy which has undergone scattering with some inhomogeneity 
in the form of a flat plane. The two most commonly considered types are sea-surface 
reverberation and bottom reverberation. Several types of volume reverberation, however, 


can be effectively modeled with an equivalent surface reverberant layer. Perhaps the most 
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well-known example of this method is the so-called deep scattering layer, which is in fact 
a thick layer of marine life which tends to vary in depth depending on the time of day 
and the season. Due to its plane-like structure, it is typically more efficient and equally 
accurate to model its effect using a surface of equivalent scattering strength located in the 
proximity of the layer volume. 

Due to the statistical nature of the large number of scatterers contributing to the 
total reverberation level at any one time, reverberation 1s usually discussed in decibels 
referenced to 1 pPascal (dB re 1pPa) instead of describing it as a time series in terms of 
pressure for a particular medium. As such, the most general form of the bistatic 


reverberation equation is: 


- IL 


scatterer -receiver 


RL = SL - TL 


source -scatterer " IS (3.1) 

SL is the source level of the projector used for the active sonar system and is typically 
given in decibels referenced to 1 pPascal at one meter (dB re 1pPa @ 1 m). For sources 
which are beamformed, this term may be different depending on the angle from which 
the energy leaves the source. SL may be given as the total level over a band of 
frequencies or as a narrowband power spectral density referenced to a 1 Hertz bandwidth 
(units of dB re 1 pPaí/Hz (2 1 m). The two TL terms represent the transmission 
(propagation) loss from the source to the scattering inhomogeneity (scatterer) and from 
the scatterer to the intended receiver of the sonar. Energy traveling from the source to 


the scatterer on one "path" will, in most cases, be scattered back to the receiver along 


many paths connecting the scatterer and the receiver. For this reason, the propagation 
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modeling which is used both in prediction and data analysis methods can be the most 
significant choice that is made in the reverberation modeling process. 

The last term, 7S,, (reverberant target strength) is the ratio (usually given in 
decibels) of incident power arriving at the scatterer on one path to the power which is 
retransmitted along a return path to the receiver. This term is dependent on both the 
characteristics of the volume or surface causing the scattering and the total area or volume 
of the scatterer. There are two possible ways in which such interactions may be 
predicted. The first, discrete backscatter prediction (Medwin, 1981, Clay et al., 1977, and 
Dyer et al., 1993), involves a detailed knowledge of the scattering surface. Each scatterer 
("wedge", "plane strip", or "facet", respectively) must be identified and its effect on 
incoming energy quantized in order to effectively predict the resulting backscattered field. 
Every interaction of incoming energy in a given propagation mechanism (eigenray or 
mode) produces an impulse in the backscattered field which is both weighted by a 
magnitude and phase term and shifted to a time which is the sum of the travel times from 
the source to the scatterer, and back to the receiver. This method is analagous to 
estimating the impulse response of the ocean as if it were a linear, time-invariant filter 
(Ziomek,1993). The total output reverberation signal is thus determined by a convolution 
of the input signal (in the time domain) with the impulse response (also in the time 
domain). As indicated by the discussion, these methods require such precise knowledge 
of the ocean environment, which is in fact time-varying, as to prohibit accurate prediction. 
The distinct advantage of this approach, however, is that it produces a signal estimate 


which is corrected for phase differences due to propagation and scattering, and hence may 


Z1 


show finer detail (higher resolution) than do the "envelope" methods discussed below. 
Methods which employ discrete backscatter prediction have been developed by Medwin 
(using Biot and Tolstoy’s "wedge" theory, 1981), Clay, et al. (using "plane strips", 1977), 
and Dyer, et al. (using "facets", 1993). 

The second method, reverberation envelope prediction, entails estimating the 
scattering strength (S, or S,) of the surface causing the reverberation depending on the 
angle or wave number of the incoming and outgoing energy and the type of surface 
and/or its estimated roughness. A great deal of research effort has been expended to 
determine scattering strength empirically for every area in the world over a wide range 
of frequencies. S, and S, are given in decibels per unit volume or area (Sy, is given in 
decibels referenced to 1 meter’ (dB re m?), S, in decibels referenced to 1 meter? (dB re 
m’) ). They are in turn multiplied by the volume or area of the scatterer (by adding 10 
log(Volume) or 10 log(Area) in decibels) to give the reverberant target strength 7S, of 
the entire scatterer. Only the amplitude level (i.e. not the phase) of the reverberant signal 
is predicted, hence the name envelope prediction. The equations which generate the 


reverberant target strength are thus: 


TSp = Sy + 10 log(Volume) 


3:2 
- $, * 10 log(área) d 


In the BSPFEX, the frequencies of concern are determined by the bandwidth of our 
receiving array, 50 to 500 Hz. In this frequency range the effects of volume reverberation 


can be assumed negligible after 1 second following the source impulse (NRL 8721, 1991). 
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For reasons discussed at length in Chapter II, the first second of data will not be 
analyzed. Hence, volume reverberation is not a concern. 

This thesis will be mainly concerned with surface reverberation caused by the sea 
bottom, since the majority of the reverberation in the BSPFEX is assumed to be due to 
bottom backscatter. This assumption is based on the oceanography in the Barents Sea and 
the frequencies which M under investigation and the Barents Sea environment. The 
effect due to the oceanography stems from the fact that the sound speed profiles in the 
Barents Sea are generally downward refracting with a strong mixed layer, both of which 
minimize sea surface reverberation and enhance bottom reverberation. 

In general, for a surface to be effective at reradiating a particular frequency of 
sound, it must have features whose size are on the order of a wavelength of sound at that 
frequency. For example, a sea surface with rms waveheight of 0.5 meter will be an 
effective scatterer for 10 KHz sound, but a weak scatterer at best for 100 Hz sound. 
Since the seas were relatively calm during the experiment (less than 0.1 meter rms wave 
height), surface backscatter can be assumed negligible compared to the effect of the 
bottom. Experimental data measured by Chapman and Harris support this assumption at 
a frequency of 500 Hz, placing the surface backscattering strength an order of magnitude 
lower than our bottom backscattering strength, with the trend serving to minimize the 
effect of surface reverberation further at lower frequencies (Urick, 1983). 

An additional effect of the bottom is subbottom AEA the reradiation or 
reflection of energy from inhomogeneities below the water/sediment interface. This effect 


can also be described as bottom volume scattering. Subbottom reverberation is not a 
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concem in the Barents Sea, however, due to the extremely thin sediment layer described 


in Chapter II and the "fast" rock bottom immediately below the sediment. 
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Figure 3.1: SUS 3 ground wave and waterborne arrivals. 


To verify this, the speed of sound in the bottom can be calculated using the SUS 
data. At the leading edge of the received signal is the ground wave, energy which has 
entered the bottom (at the critical angle for the least time path), propagated through the 
bottom to the array, and then reentered the water column (again at the critical angle) to 
be received at the array. The ground wave and waterborne arrivals for one of the SUS 


charges is shown in Figure 3.1. The process by which the difference in travel times is 
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used is displayed graphically in Figure 3.2. The speed of sound was found to be greater 


than 3000 m/sec in the bottom, verifying that the subbottom scattering is insignificant. 
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Figure 3.2: Use of ground wave to determine sound speed in the bottom. 

An additional effect is shown in Figure 3.1 immediately before the waterborne 
energy arrives. The small positive-going spike approximately 20 msec. prior to the direct 
path blast arrival is most likely energy which was trapped in the thin sediment layer. 
Using the travel time difference between it and the waterborne energy, the sound speed 
of the sediment layer is calculated as C.odimen = 1668 m/s. This number can then be used 
along with the frequency content of the trapped energy to determine the sediment 
thickness. The dominant frequency present in the energy trapped within the sediment 


“waveguide” is 286 Hz (accounting for the energy still propagating in the rock basement 
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ground wave). This yields a sediment thickness measurement of 5 meters, which is 
consistent with tabulated data for the region. The significance of the sediment sound 
speed and thickness lay in the fact that the energy trapped within the sediment layer is 
minimal, and the sound speed of the sediment is such that the layer is almost acoustically 
transparent. 

Hence, with sea-surface reverberation and volume reverberation (both in the water 
column and in the bottom) insignificant with respect to bottom reverberation, a discussion 


of bottom backscattering theory is warranted. 


B. BOTTOM SCATTERING THEORY 

The bottom is an effective reflector and scatterer of sound and acts to redistribute 
a portion of the sound incident upon it back into the ocean above it (Urick, 1983). This 
redistribution of energy is illustrated in Figure 3.3, which shows the effective "beam 
pattern” of the sea bottom for reradiating incident sound energy. Four cases are shown, 
illustrating the four general cases which occur in bottom scattering. The first and second 
general cases (cases (a) and (b)) describe bottoms which have a high impedance contrast 
with the water volume, where the acoustic impedance of a medium is defined as the 
product of the medium's density p and sound speed c. Any sufficiently hard bottom fits 
this definition, be it rocky or hard sand. Note that almost all of the incident energy is 
redistributed in the water volume above, with roughness determining the relative strength 
of the specular reflection. Cases (c) and (d) show the scattered field for bottoms which 


have low impedance contrast with the water column. Bottoms of this type are sometimes 
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Figure 3.3: Directional patterns of the scattered sound from the bottom for 
different conditions of roughness and impedance (Z=pc) contrast. (Urick, 1983) 


called acoustically transparent in that most energy is propagated into the bottom, leaving 
relatively little to be distributed back into the water volume above. Examples of these 
bottoms are soft mud or silty bottoms with a thick sediment layer. Note that bottom 
roughness is still a dominant factor in the magnitude in both the specular and scattered 


directions. 
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This brings up an interesting fact. Although the type of bottom sediment deposits 
may serve as a first-cut means of classifying bottoms in terms of acoustic backscattering, 
the roughness of the sea bottom appears to be the dominant factor in determining the 
backscattering characteristics. Just as wind speed has for many years been used as an 
indicator of surface roughness and bubble density near the surface, so is the bottom type 
merely an indicator of bottom roughness. (Urick, 1983) 

Many attempts to quantize and mathematically model bottom backscattering are 
documented. Of these, three theories of bottom scattering will be considered (as 
discussed by McCammon, 1993): 

1. Composite Roughness. 
2. Lambert’s Law. 


3. Omnidirectional Scattering. 


1. Composite Roughness 
The first theory is that of composite roughness, which McCammon describes 
as, "the most rigorous theoretical approach" (McCammon, 1993), since it can be derived 
theoretically using the Helmholtz integral equation. Composite roughness theory is 
composed of two types of scattering mechanisms: facet scatter and Bragg scatter. 
a. Facet Scatter 
For backscatter due to incident angles greater than 70° grazing, and 


forward scatter due to all incident angles, facet scatter theory is applicable (McCammon, 


1993, and Tolstoy, et al., 1987). It assumes that the bottom is composed of "facets" or 
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small plane-like structures that only specularly scatter coherent and incoherent sound 
energy. Hence, in order for a facet to scatter energy from a given raypath into a scattered 
path, the facet must be tilted sufficiently to provide a specular reflection. These types of 
models require information regarding the bottom’s rms roughness height and the spatial 
correlation length of that roughness. Most assume the surface’s roughness is Gaussian. 
McCammon presents a formula for calculating the scattering strength based on the surface 


slope and rms slope: 


Sic. =10 logl4e Cut) (3.3) 


fac 


0, 1s the local surface slope, 0,,, is the rms slope and is given by: 


VH, (3.4) 
wp 





0 


where H,,, is the rms surface roughness and L is the spatial correlation length. A is a 
scaling factor which varies depending on the surface composition. This model has been 
successfully applied in the modeling of forward scatter and has been included in the 
BISSM2 Modified Bistatic Scattering Strength Model. (McCammon, 1993) 

This model appears quite similar to the low-frequency scattering function 
for a rough surface with an isotropic Gaussian correlation function developed by Tolstoy 
and Clay (1987). This function only considers the incoherent scattered radiation relative 


to the incident radiation: 
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Sy = Dog (=) |RP. y^o? f*(0) K? (rj. e t en | (3.5) 


R is the Rayleigh reflection coefficient of the surface, and for angles less than critical, 
has magnitude equal to 1. y is the vertical wave number difference for the incident and 


scattered rays and is given by 
= —K(cos0, +c0s0 3.6) 
SS ¡ +c0s0,) 


where 0, and 0, are the incident and scattered angles relative to the vertical and k is the 


wave number given by 
i - 2f (3.7) 
C 


depending on frequency f in Hertz and speed of sound c in meters per second. © is the 
rms roughness of the surface in meters. f(O) is also an angle dependent term which, 


assuming a two dimensional problem, becomes: 


| * cos0,cosO, -sinðsinð, 





f(0)= (3.8) 
cos0, +cos8, 
r, is approximately the correlation length r, of the bottom, such that 
] 1 2 
um um em (3:9) 
ror KR 


where R range of the receiver. In the far field, r,' = r,. Finally x is the horizontal wave 


number difference, given in the following equation: 
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x = k(sinO, - sin0;) (3.10) 


Substituting these definitions into the original equation, we see the somewhat more 


familiar form: 


Sy= 10 lod (5) kto? [1+cos(8,+8,)f Fo eE | (3.11) 


Compared to McCammon's formula for facet scattering strength, we see an exponential 
term which depends on the incident angle of radiation and the spatial correlation length 
of the bottom. The A parameter in McCammon's formula appears to be equal to the 
remainder of the equation. These two equations (3.3 and 3.11) are significant because 
they were to be used to check whether the backscattering strength model used in the 


prediction routine is consistent with the approximate bottom roughness in the Barents Sea. 


b. Bragg Scatter 

Bragg scatter is applicable to low grazing angles, where composite 
roughness computes scattering based on interactions with resonant components in the 
roughness of the surface. Hence, for those components (facets) to be significant, they 
must have spatial wavelengths equal to twice that of the incident sound. This concept of 
resonances is similar to that seen in diffraction gratings, hence the name Bragg scatter. 
The input parameters for this model are the surface height power spectrum W(K) at the 
resonance wavenumber K = Kx (Eq. 3.10), which can be derived from the surface rms 
roughness and correlation length, given the surface roughness has a Gaussian distribution. 


Bragg backscattering strength in two dimensions is given by 
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S rase" vog SOD fey eat UD (3.12) 
x 
where the angles 0, and 0, are the grazing angles of the incident and scattered energy, as 
in Eq. 3.3 (McCammon, 1993). Again, this theory depends on knowledge of the assumed 
wavenumber dependence of the spectrum W, which in general is not well documented for 
most of the ocean bottoms on the planet. This model is based solely on a statistical 
model of the physical structure of the reflecting surface. Hence, effects due to 
deterministic structures (bottom features), volume scattering effects (subbottom scattering), 
or refraction (the basis of discrete backscatter models) are not accounted for. As a result, 


the effective fourth power dependence on the sine of the grazing angle causes S,,,,, to be 


inaccurate (low) below 20°. 


2. | Lambert's Law 
The second theory is Lambert's law of diffuse scattering. As noted by Urick, 
Lambert's law is, "a type of angular variation which many rough surfaces appear to 
satisfy for the scattering of both sound and light" (Urick, 1983) Noting the original 
application of this "law", he continues: 
Although many materials follow Lambert's law closely in scattering light, none 
does so exactly. Lambert's law applies specifically to the radiation of light by 
radiant, absorptive materials; the "law" should properly be called Lambert's "rule" 
for scattering. Nevertheless, it is, as we have seen, a good description of the 
backscattering of sound by very rough bottoms. (Urick, 1983) 
The original application of Lambert’s law was to optics, in which it is known 


as Lambert’s cosine law of diffuse reflection. As described by Houstoun (1938), a 


portion dS of a surface has apparent bnghtness which vanes with the sine of the grazing 
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Figure 3.4: Geometry for discussion of Lambert’s law of diffuse reflection. 


angle of the incident light, but which does not vary with the grazing angle of the reflected 


light. Specifically, 
I _ B Lincidens Sin(8)) (3.13) 
apparent r2 


where, as shown in Figure 3.4, p is a maternal dependent constant, J;,-isen is the intensity 
of the incident light, 0, is the grazing angle of the incident light and r is the distance 


between the surface and the observer. 7,,4,,^ I1 dS sin(0j) is the total light received by 
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the area dS due to light of intensity J. A practical example of this relationship is a 
luminous sphere, for example, a red-hot copper sphere heated over a flame. When 
viewed from a distance, it appears as a disc of uniform brightness over the entire surface 
when viewed from any angle. Hence, the radiation from an element of the surface 
making an angle 0, with that element of the surface must consequently be proportional 
to sin(8,) (Houstoun, 1938). The name "cosine law" comes from the fact that in its 
original development, angle were measured with the normal to the surface, not the grazing 
angle with the surface itself. 

In the application of Lambert’s law to acoustics, the scattering strength varies 
as the product of the sine of the incident and -scattered grazing angles times a constant u 
that is a function of the material properties of the scattering surface. Thus, scattering 


strength is given by: 
S Lambert= 10 log(psin6,sin6,) = pgp +10 log(sin6,sin8,) (3.14) 


If all of the incident energy were redistributed in the water column, with none lost by 
transmission into the bottom, it can be shown by integration that p = I/n, or 10log(pu) = 
Hap = -4.97 dB. In 1961, K. V. Mackenzie reported success using Lambert’s law with Hag 
= -29 dB to model bottom backscattering in deep water at 530 and 1030 Hz. His 
conclusions are not viewed as a rigorous theoretical result, but instead a simple 
geometrical argument which fits the angle dependence of "a lot” of measured bottom 
backscattering data at moderate grazing angles, specifically 20% to 60%. (Mackenzie, 


1961) 
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Figure 3.5: Reverberation following a 2 1b. explosive charge detonating at 244 
m in water 1980 m deep, as observed with a nearby hydrophone at a depth of 41 
m. Filter band 1 to 2 kHz. (Urick, 1983) 


The choice of 20° to 60° as limits for the application is not arbitrary. Instead, 
they result from the geometry of the experiments that were conducted in deep water and 
the method by which the results were analyzed. As shown in Figure 3.5, reverberation 
due to the first surface reflections, volume scattering, scattering in the deep scattering 
layer, and bottom backscatter are separated in time. At approximately 2.4 seconds, the 
first bottom returns, the so-called fathometer returns, come in. At 2.7 seconds the 
fathometer return from the surface reflection of the explosion arrives at the receiver, 
followed by a long (~2 second) bottom reverberation return. At 2.7 seconds, the 
approximate grazing angle of the energy interacting with the bottom is 60°. At 4.8 
seconds, after which the next set of bottom-surface reflections arrive at the receiver, the 
bottom grazing angle of the direct path energy (energy which has traveled through the 


water column only, with no interactions with the surface or bottom) has decreased to 
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approximately 20°. Hence, the limits for the applicability of this monostatic reverberation 
measurement arise naturally from the experimental procedure. The term monostatic is 
appropriate in this context because the only scattered paths which are significant are those 
which follow the same grazing angle as the incident energy. 

The bistatic dependence of Lambert's law has been examined by P. B. 
Schmidt, again in deep water, in a 1971 paper. He found good agreement with bottom 
scattering data for incident angles from 3.5? to 78? with scattered return angles from 21? 
to 84° (McCammon, 1993). The value of p4g must be set empirically, however, which 
requires considerable analysis in finding the variation of Hag With frequency, azimuth 
angle, and sediment roughness. This theory has been successfully implemented by J. W. 
Caruthers, et al. in the Modified Bistatic Scattering Strength Model (BISSM2) 
(McCammon, 1993). They found p4,= -31.5 dB to best fit their data. 

In shallow water, the problem becomes somewhat more complicated. Figure 
3.6 shows the return from a 1.8 lb charge used in the BSPFEX. Instead of having clearly 
separated volume, surface, and bottom returns, they instead all come in simultaneously. 
The implications of this are two-fold. First, the dominant scattering mechanism must be 
used to predict the reverberation level (in this case, bottom reverberation). Second, the 
problem is truly bistatic in nature, since the travel time difference between direct path 
rays and rays which have multiple interactions with the surface and bottom are small, on 
the order of the decay time of the bubble pulse of the explosion (see Section C for a 
discussion of the bubble pulse and other underwater explosion phenomena). As a result, 


for each bottom interaction, there are multiple scattered angles which significantly 
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Figure 3.6: Reverberation following a 1.8 lb explosive charge (Mk 82 SUS) 


detonating at 91.4 m in water 275 m deep, as observed with a hydrophone at a 
depth of 121 m. Filter band 50 to 500 Hz. 


contribute to the reverberation. Only one of these multiple scattered angles is the 
monostatic return; the scattering is thus bistatic in elevation. As will be shown in Chapter 
IV, the scattering is in fact bistatic in azimuth also. 

In Chapter IV, the method used to predict bottom reverberation in the Barents 
Sea will be based on an implementation of Lambert’s law in shallow water taking 
advantage of the vertically bistatic nature of S,,,,,,. This problem could not be attacked 


in this manner in the past because of the lack of adequate computing power. 
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3. Omnidirectional Scattering 
The general form of the omnidirectional scattering theory, also known as the 


Lommel-Seeliger law, is (McCammon, 1993): 
Somni = 10 log( psin®,) = p 5 * 10log(sin0,) (3.15) 


Again, Hag is dependent upon the bottom sediment properties, frequency, and azimuth (of 
the scattered energy). This theory differs from Lambert’s law in that it only depends on 
the sine of the incident energy’s grazing angle. McCammon (1993) argues that at low 
grazing angles (less than 30°), both monostatic and bistatic data support the first power 


dependence on sin(0,). This theory will be tested in Chapter IV. 


C. UNDERWATER EXPLOSION THEORY 
Underwater explosions have been used for many years as acoustic source in many 
diverse applications (Weston, 1960): 


* Geophysical prospecting, especially by the oil industry (where their use has been 
quite successful). 


e Underwater signaling, including distress signaling. 


e Sound sources for explosive echo ranging against submarine targets (though such 
use is not popular among submariners, including the author). 


* Underwater sound propagation experiments, including reverberation studies. 
In the latter application, they have become the de facto standard source due to their 
energy density, near impulse-like characteristics, and low cost. They also have distinct 
advantages over other impulse-like sources (such as sparkers) in that they are easy to 


deploy and are not limited in depth. This is not to say they are perfect, however, since 
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the signal they produce is somewhat random and is riddled with distinguishing 
characteristics which vary depending on the type of explosive/detonator, the fusing 
mechanism, the size of the charge, and the depth of detonation. Since they have been 
used for so long, many of these characterisucs have been measured and quantzed, and 


hence they have become relatively predictable sources of sound. 


1. The Underwater Explosion and Gas Globe 

When a high explosive is detonated, the detonation wave propagates through 
the mass of explosive material, converting it to an incandescent gas at a very high 
pressure. In high explosives like TNT, this detonation wave or "shock" wave is 
supersonic with velocity between 5,000 and 10,000 m/s. A spherically symmetric shock 
wave is then radiated into the surrounding water with a pressure rise that is practically 
instantaneous, followed by an exponential decay with a time constant on the order of a 
fraction of a millisecond. The gas bubble or "globe" created by the explosion grows with 
the shock wave, but due to inertia it overshoots its equilibrium radius, and hence the 
radiated pressure becomes slightly negative as the bubble begins to contract. As it 
contracts it again overshoots its equilibrium radius, and when it reaches its minimum 
volume, a positive pulse with strength and duration comparable to the initial shock wave 
is developed. Successive pulses grow increasingly weaker. This damped radial 
oscillation continues, creating positive pulses at each point of minimum radius, until the 
globe either reaches the surface or breaks up into smaller parts. The oscillatory signal 
of these explosions is known as the bubble pulse. The pressure signature of an 


underwater explosion at short ranges consists of the shock wave followed by a small 
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number (less than 10) of bubble pulses. At longer ranges or in shallow water, the 
received signal is complicated by multipath effects which create more positive and 
negative impulses for each path, with surface reflections causing negative-going impulses. 


(Weston, 1960, and Urick, 1983) 


2. The Bubble Pulse 

The time difference between the shock wave and the first bubble pulse is 
determined by the weight of the charge, the type of explosive, and the depth at which it 
is detonated. This time difference T is predictable and can be found theoretically (Urick, 
1983): 

. Kw" (3.16) 
(d+33)% 

K depends on the type of explosive (K = 4.36 for TNT), w is the charge weight in 
pounds, and dis the detonation depth in feet. The reciprocal of T is known as the bubble 
pulse frequency and is measurable from data. For a Mk 82 SUS charge detonating at 300 
ft, the predicted bubble pulse frequency is 23 Hz, which agrees well with the Barents Sea 
data. This measurement is important, since it proves that the SUS charges were 
detonating at the correct depth with the correct strength in the shock wave. 

Energy is radiated as long as the bubble continues to oscillate, but for large 
charges and shallow depths (i.e. when the bubble is large), the bubble rises appreciably 
during one oscillation, resulting in a new source depth and a different multipath problem. 
Additionally, it affects the energy spectrum of the source, since the frequency of 


oscillation will shift as the bubble migrates toward the surface. Because of the possible 
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implications of this bubble migration, it is important to know when the effect is 
significant. Weston (1960) proposes that the critical depth below which this effect is 


small is given by 


d, - 200W/? (3.17) 


where W is the charge weight in pounds. This corresponds to 232 ft for a 1.8 lb charge. 
Since the Mk 82 SUS detonated at 300 ft, bubble migration is not a concern in the 


BSPFEX. 


3. The Energy/Power Spectral Density 

Since the pressure signature of a SUS is by its very nature random, it is 
necessary to move to the frequency domain in order to avoid complications caused by 
modeling such an explosion in the time domain. This shift also entails a choice of either 
energy flux density or power spectral density of the shock wave and first few bubble 
pulses (which contain almost all of the energy generated by the explosion). This thesis 
will use power spectral densities in units of decibels referenced to 1 pnPa^/Hz (2 1 meter. 
Weston and Unck both provide spectra for charges. Weston’s work, while excellent, is 
somewhat dated. Urick includes a disclaimer with his spectrum calculation in which 
describes his curves are, "based upon, and are useful only for, broad-band levels, and are 
misleading for narrow-band work in not showing the fine-frequency spectral variation 
caused by the multipulse character of an explosive signal.” (Urick, 1971) Hence, a more 


recent source 1s needed. 
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Comparison of Mk 82 SUS Source Level Values 
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Figure 3.7: Comparison of Mk 82 SUS source level spectra. 


Subsequent to the above researchers, source level measurements for various 
types of charges were made by Stockhausen (1964), Christian and Blaik (1965), Turner 
and Scrimger (1970), Gaspin and Shuler (1971), and Chapman (1987). Prior to Chapman, 
the most widely used values for the source levels of shallow SUS were Gaspin and 
Shuler’s. Their estimates were derived from hand-drawn waveforms that were based on 
experimental measurements of the waveforms of deep explosions extrapolated to shallow 
depths. Between Gaspin and Shuler’s work and Chapman’s, there were several 


measurements performed that were in general agreement with the Gaspin and Shuler 
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predictions, but there were considerable differences among the data sets. (Chapman, 
1987) 

Chapman’s results are presented in 1/3 octave bands from 5 to 630 Hz. At 
least 30 charges (1.8 lb SUS) were averaged for each detonation depth reported. The 
standard deviations were within 1 dB for most frequency bands (from 50 to 500 Hz all 
standard deviations were less than 1 dB). Chapman notes that actual explosion depths 
were consistently deep in all measurements. He theorized, with support from a previous 
experiment, that this was due to the relatively cold water (9° C) in which the experiment 
was conducted. Instead of the charges averaging the nominal 91.44 m (300 ft), he noted 
an average depth of 99.6 m (326 ft). This effect should be insignificant in the BSPFEX, 
as the effect on bubble pulse frequency is slight. Also, since the propagation model 
assumes a monostatic problem, a deeper detonation actually moves the detonation closer 
to the top of the array, making the monostatic assumption more viable. 

Chapman’s results agree well with Gaspin and Shuler’s, with the possible 
exception of the lower frequencies. Chapman points out clearly, however, that this effect 
was expected because Gaspin and Shuler did not specifically account for the impulse of 
the waveform, which would throw off their low frequency predictions (Chapman, 1987). 
In this thesis, the source levels measured by Chapman will be used, utilizing interpolated 
values (cubic spline fit) to obtain an approximate source level every 12.5 Hz from 50 to 


500 Hz. This is shown by the dashed line in Figure 3.7. 
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In the following chapter, the SUS model will be used as part of the input for 
a reverberation envelope prediction system designed for the Barents Sea. The predictions 


made by the system will then be compared to the shot data to validate its performance. 


IV. METHOD AND RESULTS 


A. PREDICTION METHOD 

As described in Chapter I, propagation modeling for the Barents Sea has been 
performed using ray traces provided by HARPO. The propagation model from HARPO 
is then used along with the SUS source model descnbed in Chapter III and either 
Lambert’s law or omnidirectional scattering bottom backscattering functions to produce 
a predicted reverberation envelope. Each of the components in the prediction routine are 


described in the sections below. 


1. HARPO 

HARPO (Hamiltonian Acoustic Ray-tracing Program for the Ocean) is a 
FORTRAN computer program that traces the three-dimensional paths of acoustic rays 
through model oceans by numerically integrating Hamilton's equations, which are a 
differential expression of Fermat’s principle. It was developed by R. Michael Jones, J. 
P. Riley, and T. M. Georges at the Wave Propagation Laboratory in Boulder, Colorado 
(an Environmental Research Laboratory of the National Oceanic and Atmospheric 
Administration). As described in Jones, et al. (1986), it uses continuous models for the 
sound-speed distribution and bathymetry in the ocean, hence avoiding the false caustics 
and discontinuous raypath properties encountered in conventional ray-tracing methods. 


Only geometrical effects are accounted for in the ray-trace; that is, no diffraction or 
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partial reflection corrections are applied by HARPO itself. Amplitude along a raypath is 
not explicitly computed; that is left to a post-processing program which accounts for 
boundary interactions and losses due to ray tube spreading and absorption. The program 
does produce a step-by-step account of each ray’s progress, including all turning points 
and boundary interactions, and outputs the data in machine-readable form for post- 
processing. The post-processing tool chosen for this thesis was MATLAB. 

The inputs for HARPO in its present form (on the Naval Postgraduate 
School’s Cray Y-MP EL98) are a data file which contains the three-dimensional sound 
speed field in the region (sound speed at each longitude/latitude/depth point), a data file 
which contains the two-dimensional bathymetry for the same region (depth at each 
longitude/latitude pair), the source/receiver location, the desired launch elevation angles, 
and the desired launch azimuth. It then computes closed-form expressions for the 
bathymetry and the sound speed field using a user-specified number of "empirical 
orthogonal functions" (EOF's) to model the sound speed field (Chiu, et al., 1993). The 
number of EOF's required depends on the original number of sound velocity profiles 
(SVP’s) entered and the maximum error desired in the continuous representation of the 
field. 

It should be noted that HARPO cannot directly compute the raypaths that 
connect a specified source and receiver (or scatterer in the case of reverberation). To find 
such eigenrays, one usually launches a fan of rays at small usce in elevation and/or 
azimuth angle and linearly interpolates for the rays that reach the desired receiver 


location. In the specific case of bottom reverberation, since each bottom interaction is 
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Latitude 





Longitude 
Figure 4.1: CTD locations and extrapolated points along isbaths. 


identified in HARPO's output "rayset," this turns out to be an interpolation along the 


bottom only (assuming the bottom is locally flat). 


The bathymetry for the Barents Sea was obtained using the WORLDBATHY 
database in the Naval Postgraduate School’s Oceanography laboratory. The data 
originally came from DBDBS, a standard U. S. Navy world bathymetric database. The 


sound speed field was measured as part of the BSPFEX. A total of 72 
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Conductivity/Temperature/Depth (CTD) measurements were taken, 56 of which were 
within 38 km of the array and hence were usable in the reverberation experiment. To 
cover the entire area, the CTD measurements were copied to new locations along isobaths, 
smoothed using an 11 point Hamming smoother, and then an inverse distance algorithm 
(MATLAB, 1992) was applied to interpolate the sound velocity field at regular grid 
locations for 10 — and 9 latitudes around the array. For the three-dimensional 
sound speed field in the Barents Sea experiment area, a total of 90 SVP's were used to 
estimate the sound speed field. To keep error less than 0.1 m/s, 26 EOF’s were required. 

To estimate the horizontally bistatic nature of the experiment (the source SUS 
and receiving array were -1 km apart) using a monostatic propagation model, a 
source/receiver position that was in the center of the SUS drop points and array location 
was used. From that point, rays were launched on several azimuths over a fan of 
elevation angles. The runs were originally performed using a three-dimensional 
environment, but it was quickly discovered that doing so violated the assumptions made 
in modeling the reverberation using an N x 2D approximation. To correct this 
discrepancy, the final HARPO runs were performed using two-dimensional "slices” from 
the three-dimensional environmental data. Four launch azimuths (000, 090, 180 and 270) 
were used, with rays traced every 0.5° in elevation from -60° to -20° and from +20° to 
+60°, and every 0.1° in elevation from -20° to +20°. The results of both the 2D and 3D 


data runs are described later in this chapter. 
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2. | Post-processing of HARPO: rbreakbr.m 

rbreakbr.m is a MATLAB script which reads in the output data file from 
HARPO and processes it to find eigenrays between the source/receiver location and 
desired range points along the azimuth on which HARPO was run. For each eigenray to 
a given bottom point, it determines the ray's amplitude, path length, horizontal range, 
depth, phase, arrival angle and travel time. The local sound speed at the bottom, original 
launch angle of the ray, and vertical distance between it and the adjacent ray used for ray 
tube spreading calculation are also stored. 

Ray path lengths are determined by simply summing the straight line distances 
between HARPO rayset output points. Stepping through all output points for each ray, 
boundary interactions and turning points are accounted for individually. At each, the type 
of interaction is identified and numbered using a system described by Franchi et al. 
(1984). Each interaction or turning point is numbered as an order contour, with all 
surface interactions/upper turning points receiving odd numbers, and bottom 
interactions/lower turning points even. If the first interaction is a bottom interaction, it 
is hence identified as #2. Each is also identified as to the type of point (boundary 
interaction or turning point). This numbering method turns out to be key in later 
processing when interpolation is performed between adjacent launch angle rays at bottom 
interactions with the same order contour. 

Amplitude calculations for a ray at a given bottom interaction are determined 
by the combined effects of surface and bottom interactions the ray has experienced up to 


that point and ray tube spreading. At each surface interaction, the loss mechanism 
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assumed is due to scattering loss only using the Eckart coherent reflection coefficient Rc 


such that 


R.-8 e 27 oon (4.1) 


where K, the Rayleigh reflection coefficient, is equal to 1, O,uyace 15 the rms waveheight 
in meters (less than 0.1 meter during the experiement; O meters used in analyses), and y 
is the horizontal wave number (Tolstoy, et al., 1987). For each bottom interaction, the 
loss mechanism assumed is due to scattering loss using the Eckart coherent reflection 
coefficient along with SX. For angles less than critical (68° for the bottom in the vicinity 
of the array, as measured using the ground wave of the SUS charges) the magnitude of 
the Rayleigh reflection coefficient is 1. Rays which have angles higher than this are so 
severely attenuated that their contribution to the reverberation level is insignificant. 
Hence, the magnitude of the specularly reflected ray after a bottom reflection is also 
determined by Equation 4.1, where 64,15. 15 the rms roughness of the bottom in meters 
(Tolstoy, et al., 1987). Obəætm Was determined from fathometer recordings to be 
approximately 0.5 meters. At the frequencies of interest in the BSPFEX, loss due to 
absorption in the bottom is negligible. Ray tube spreading is calculated using 


conservation of power within the ray tube given by the following equation: 
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cos(8,) (4.2) 
cos(0) 


P E A0, pc 
5 Th PoCo 








P is the acoustic pressure at horizontal range r from the source, P, is the acoustic pressure 
l meter from the source, AQ, is the angular difference, in radians, between adjacent 
launch angle rays, h is the vertical distance between the two rays at range r, p and p, are 
the water densities at the bottom point and at the source (assumed equal), c and c, are the 
sound speeds at the bottom point and at the source, 9, is the launch elevation and @ is the 


elevation angle of the ray at the bottom interaction (Clay, et al., 1977). 


Method of Interpolation Between Adjacent Rays 
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Figure 4.2: Notional view of interpolation between adjacent launch angle rays. 
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Once all bottom interactions have been identified and eigenray data for each 
has been stored, the program then steps out in range to the points specified (at regular 
intervals) and finds adjacent launch angle rays which have bottom interactions with the 
same order contour which "straddle" the desired bottom point. This is illustrated 
graphically in Figure 4.2, where the two solid lines represent adjacent launch angle rays 
(which have bottom a with the same order contour) which straddle the location 
of a desired output point on the bottom. A linear interpolation is then performed in 
range, assuming the bottom to be locally flat, to find ray amplitude, path length, depth, 


phase, travel time, arrival angle and local sound speed at the desired point. 


3. Combining Eigenrays, SUS Model and Scattering Function: rlmaker.m 

The goal of the rlmaker.m routine is to estimate the reverberation level 

prediction algorithm suggested by Franchi, et al. (1984) and more recently by Bucker, et 
al. (1993). They both suggest a double integral over the reverberating area with multiple 
elgenrays connecting each "patch" dA on the bottom with the source/receiver. In order 
to model reverberation in this manner using a true three-dimensional model, ray traces 
would have to be done from each SUS location over a fan of launch elevations for each 
azimuth in a 360° fan surrounding the SUS. A similar run would have to be done for 
each of the other three SUS and for the vertical line array. Then, rays connecting the 
SUS and the vertical line array with a given patch on the bottom would have to be picked 
out of the set, and the appropriate amplitude/travel times -—5 to each. This is a 
formidable task, yet Bucker, et al., claim to have successfully modeled reverberation in 


this way with some simplifications in their RUMBLE product (Bucker, et al., 1993). 
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Confining the problem to a two-dimensional slice of the ocean and assuming 
a horizontally monostatic problem lessens the computational load significantly, and 
enables the routine to run within a reasonable amount of time. It also ensures that the 
rays connecting the source to a given scattering patch are the same as the rays which 
connect that patch back to the receiver. The reverberation level at the receiver due to a 


source with intensity I(t) is thus 


B(0,5,,9 BOQ, 
Ro=ff Y Y It-T(r)-T()) oo 9? 


0(0;,0.) dA 
A iem(r) jem(r) Lr )LAr ) oat 


(4.3) 
The two T terms represent the travel times along the ith raypath from the source to the 
bottom point and for the jth return path in the set of raypaths represented by the set m(r) 
connecting the source/receiver and the bottom patch dA which is at range r. The two B 
terms represent the beamforming factors to be taken into account on the source and 
receivers and the two L terms account for losses in amplitude. O is the scattering strength 
of the bottom, depending on the incident energy grazing angle 0, and the scattered 
energy’s launch angle 8,;. No azimuthal dependence is modeled. This method is 
displayed graphically in Figure 4.3. 

Franchi et al. note that it was not possible for them to evaluate Eq. 4.3 
directly due to a lack of a three-dimensional propagation model. Bucker et al. were able 
to evaluate the intergral because they had both the model and the required computing 
assets. At the Naval Postgraduate School, there were adequate computing facilities and 


a three-dimensional propagation model, but the run time required to properly model the 
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patch on the bottom, 
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eigenrays returning 

to the array (2 shown). 








Figure 4.3: Illustration of multipath reverberation problem. 


reverberation in shallow water proved prohibitive at the time of this writing. 
Given the asssumptions stated above, and estimating the area integral by 


summing the contributions to reverberation by rings on the bottom, Eq. 4.3 simplifies to 


RL A= Y p B Kt-T,-T, f) LOL es f xrL-ri) (44 


B is now the effective beamwidth of the point source/receiver given the region for which 
the ray trace is applicable. For the four azimuth approximation to be made, B = 0.25. 
I(t,f) is the intensity of the SUS at frequency f and is assumed constant at its peak level 
for 0.1 sec. The summation over k performs a power sum of the contributions from each 
ring out to a range of ~38km, placing the reverberant energy into the appropriate 0.1 
second time bin depending on the two-way, vertically bistatic travel time of the ij raypath. 


The two loss terms and the bottom backscattering strength are frequency dependent. r, 
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and r,, represent the inner and outer radii of the kth ring which has its center at 
(r*r,)/2, and such that rings do not overlap. Hence, RL(t,f) is the predicted 


reverberation peak envelope at the desired frequency. 


B. BARENTS SEA DATA ANALYSIS 

To most closely estimate a monostatic problem, only the phone closest to the 
detonation depth, phone O, the top phone, was analyzed for reverberation level. 
Broadband peak level measurements were performed using a peak picking algorithm and 


a straight conversion from voltage to sound pressure level in decibels: 


RL(t)- 20 log( voltage x10?) (4.5) 
ollas levels were corrected for phone sensitivity problems and non-zero mean as 
described in Chapter II. 

Narrowband measurements were performed using a Short-Time Fourier Transform 
(STFT) power spectral density estimate using an N = 128 point rectangular window, with 
1/2 window overlap in segment intervals (Therrien, 1992). If the fast Founer transform 
(fft) of data vector x(t) is given by X(f), then the narrowband level of x(t) in a 1 Hz 


bandwidth is given by 


RL(f) e 2E (4.6) 


f, 1s the sampling frequency (1600 Hz). The factor of 2 is required to maintain 


conservation of power when only considering positive frequency. Prediction results and 


po 


BSPFEX data are presented in Section 4.C. Statistical analysis of the reverberant returns 


are described in Section 4.D. 


C. PREDICTION RESULTS 

Narrowband predictions were made for 1 Hz bandwidths centered at 50, 100, 200 
and 400 Hz. In the process of developing the results, the backscattering coefficient of the 
Barents Sea has been determined. Both Lambert’s law and omnidirectional scattering 
were considered. In general, the choice of the backscattering function tends to change 
the slope of the estimate with time, while the choice of the backscattering strength shifts 
the estimate up and down without slope change. Lambert’s law provides a better fit to 
the measured data. The values of pin Table 4.1 were found to best fit the Barents Sea 
data. They are consistent with the findings of McCammon (McCammon, 1993). 


TABLE 4.1: BACKSCATTERING STRENGTH PARAMETER y FOR 
SEVERAL FREQUENCIES IN THE BARENTS SEA. 












Determination of which scattering function is better is somewhat subjective, in that 
the choice of parameters in the propagation model determine which scattering function 
fits better. Figure 4.4 demonstrates this by showing the difference obtained using 
Lambert’s law and omnidirectional scattering. Any loss term which increases with range 


also changes the slope of the estimate, more dramatically, in fact, than the choice of the 


scattering function. Specifically, volume absorption, surface scattering loss, bottom 
scattering loss, and bottom absorpnon are all loss terms which tend to act in a roughly 
linear manner to increase transmission loss (TL in decibels) with range. Therefore, 
perfect, or near perfect, propagation modeling is required in order to determine which 
scattering function is more accurate. In shallow water, this may not be a reasonable 


request, especially given the nature of the multipath problem. 


Comparison of Lambert’s Law and Omnidirectional Scattering 
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Figure 4.4: Comparison of Lambert's law with omnidirectional scattering at 200 
Hz, overlaid on 200 Hz data from the BSPFEX. 
Figures 4.5 through 4.8 show the predictions and data collected during the BSPFEX 


that generated the numbers in Table 1. At 50 and 100 Hz, the propagation model is 


suspect for two reasons. The first is the lack of an adequate loss mechanism to fully 
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explain the slope of the BSPFEX data. This is especially visible at 50 Hz, as the 
prediction begins to increase at the 20 second mark. The second reason is the more 
fundamental question of whether ray theory is correct at lower frequencies. For these 
reasons, the predictions at 200 and 400 Hz are thought to be much more accurate than 
those at the lower frequencies. For the validity of these models to truly be determined 
in a general sense, a closer examination is required using data from several different 
regions. 

The trend indicated in the rising p with frequency is significant. It is also 
consistent with rough surface scattering theory as described in Chapter HI. Referring 
back to Equation 3.5, it is clear that the scattering function varies with f in a complicated 
fashion. From the predictions in Figures 4.5 through 4.8 and the values listed in Table 
4.1, there appears to be between 3 and 8 decibels per octave difference in scattering 
strength over the frequency band considered. Again, this conclusion should be tempered 


with the knowledge that the propagation modeling may not be exact below 200 Hz. 
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Figure 4.5: Reverberation level prediction and Barents Sea results at 50 Hz. 
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Figure 4.6: Reverberation level prediction and Barents Sea results at 100 Hz. 
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Narrowband Reverberation Level Predicted for Barents Sea 
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Figure 4.7: Reverberation level prediction and Barents Sea results at 200 Hz. 
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Figure 4.8: Reverberation level prediction and Barents Sea results at 400 Hz. 
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Figure 4.9: Plot showing horizontal refraction of three rays launched from the 
vertical array on azimuth OOOT. 


A significant effect not taken into account in the two-dimensional prediction 
routines is that of honzontal refraction. The significance of this effect is clearly 
illustrated in Figure 4.9. In it are shown the raypaths for three different launch angle rays 
(+25°, +40°, and +60°). The bathymetry proves to be the dominant effect on the rays as 
they refract eastward away from the shallow region to the northwest of the array. This 
Causes a two-fold problem in the reverberation model. First, the assumption that a ray 
can be scattered to a return path with some angle other than its incident angle is invalid, 


unless a fan of launch azimuths is used along with a three-dimensional propagation 


61 


model. For example, assuming that a single launch azimuth will suffice for, say, a 90° 
sector is incorrect, since this would imply that energy from a 25? launch elevation striking 
the bottom at (74.65? N, 023.66? E) would somehow be scattered into a 45? received 
elevation eigenray which emanated from (74.60? N, 024.00? E). This is simply not a 
good assumption. Of course, modeling propagation in a three-dimensional environment 
with a two-dimensional model is not a good assumption either, but at least doing so is 
consistent with the method of reverberation prediction. The second problem is the effect 
of horizontal refraction on signal amplitude due to ray tube spreading in the horizontal. 
From Figure 4.9, this problem clearly appears to be significant. Hence, to accurately 


model the reverberation in a shallow region, a three-dimensional model should be used. 


D. STATISTICAL ANALYSIS 

A statistical analysis of the BSPFEX data has been performed. Since the 
reverberant return from an explosion in shallow water is by its very nature a non- 
stationary process in time, the wideband returns (sampled at 1600 Hz) were broken into 
5000 point segments and analyzed individually, assuming approximate stationarity over 
the segment. The standard deviation 6, of each segment was determined using 
MATLAB, and the data segment was then divided by 6, to normalize the segment to 
standard deviations of pressure, a unitless quantity. This allowed segments to be 
compared directly to determine the underlying probability density function of the entire 
reverberation signal. Dividing the pressure by the standard deviation is essentially the 


Same as dividing by the rms pressure, or normalizing the signal power to one unit of 
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Normalized Probability Density Function of Reverberation 
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Figure 4.10: Normalized probability density function of reverberation signal. 
Dots show mean experimental values; vertical lines show 95% confidence 
intervals. Gaussian density overlaid for reference. 


power. The probability density function f,(x’) of each segment was then estimated using 
the histogram subroutine in MATLAB (1992) and the equation: 
Ko 
x? 5 A (4.7) 
fx ur 
where K is the number of data points used for the estimate, Ax is the size of the analysis 
bin in the histogram, and K,. is the number of data points which fell in the bin centered 


on x? (Therrien, 1992). The average over time was also complemented with an ensemble 
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average over the four SUS returns to yield the mean probability density function shown 
in Figure 4.10. The dots on the plot represent the mean probability density function 
estimate over all four SUS returns, and the vertical line through each dot shows the 95% 


confidence interval determined using 


/ 
ro? m». u, < X+ Ox nan) (4.8) 


VN VN 


where x is the sample mean of the random variable x, g’, is the sample standard deviation 
Of xX, t,.y2 15 the value of the Student ¢ distribution for n = N - 1 (N is the number of data 
points in the sample x), and 2 7 - P fora( P*100 ) percent confidence interval. 
(Bendat, et al., 1971) 

Clearly the normalized wideband SUS returns have a stationary Gaussian 
distribution (i.e. a normal density function). Individual unnormalized segments are still 
Gaussian processes, but are not stationary since they have variance proportional to the 
instantaneous power of the return signal. The erratic estimate in the vicinity of zero 
standard deviations is due to quantization noise of the A/D converter in the measurement 
system. This noise becomes significant in the normalized data late in a given SUS return 
when the data is just barely above background noise. 

In order to examine the statistics of narrowband data, the data were again segmented 
in order to compute fast Fourier transforms (fft's) of each 128 point segment. In order 
to prove that each frequency bin has a bivariate Gaussian distribution, the real and 


imaginary portions must be shown to both have Gaussian distributions. This was done 
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Figure 4.11: Narrowband probability density function estimates at 50, 100, 200, 
and 400 Hz. Dots show mean experimental values; vertical lines are 9590 
confidence intervals. Gaussian density function overlaid for reference. 
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using a method similar to that done on the wideband data described above. The results 
of this analysis are shown in Figure 4.11. As before, the dots represent the sample mean 
of the probability density function estimate over all four SUS charge returns. The vertical 
lines are the 95% confidence intervals. The narrowband data again clearly has a Gaussian 
distribution at all frequencies. 

These results are significant, because verifying the signals to be Gaussian random 
processes allows a great deal of existing performance prediction algorithms to be applied 
in the analysis of active sonar systems in the Barents Sea. The results are not surprising, 
since the data comes from a multipath acoustic problem which is by nature random due 
to the effects of the oceanographic and bathymetric features over the entire azimuth. 
Nonetheless, they are quite valuable to a potential researcher or active sonar system model 
for the Barents Sea. This, however, does not necessarily imply that these results are true 
in general. In order to extend this assumption to other shallow regions, a similar study 


must be done for that region. 
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V. CONCLUSIONS 


The objectives of this thesis were to consider several methods of reverberation 
prediction and compare these methods to the results from the BSPFEX, to determine the 
reverberant characteristics of the Barents Sea, and determine the statistical characteristics 
of the reverberation measured in the BSPFEX. The reverberation prediction routines 
which were developed for this shallow environment utilized multipath propagation 
modeling along with a vertically bistatic scattering function for the bottom. They have 
been used to accurately model the reverberation due to SUS charges dropped during the 
BSPFEX. It has been determined that Lambert’s cosine law accurately describes the 
bottom backscattering function using a parameter which varies depending on bottom 
composition and roughness. Furthermore, that parameter has been determined at several 
frequencies for the Barents Sea. Omnidirectional scattering theory, however, has not been 
ruled out as a viable reverberation prediction tool. Finally, the statistics underlying the 
reverberant returns in the shallow water of the Barents Sea have been proven to be 
Gaussian distributions, both for broadband and narrowband data. 

These results do not, by themselves, provide all the information required to predict 
the performance of a LFAA system in shallow water. They do however, provide key 
information which is necessary to achieve this goal. Knowledge of which scattering 


function to use and the parameters which drive the scattering is essential to predict 
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whether reverberation is a limiting factor. If it is, knowledge of the statistical properties 
of the reverberant returns is essential in such assessments. 

Further research effort should be directed at modeling the reverberation in different 
shallow regions using different types of sources and receivers. Such modeling should 


consider three-dimensional effects both on propagation and scattering. 
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